Finite-length soliton solutions of the local homogeneous 
nonlinear Schrodinger equation 

E.G. Caparellit V.V. Dodonovtand S.S. Mizrahi* 
Departamento de Fisica, CCT, Universidade Federal de Sao Carlos, 
Via Washington Luis km 235, C.P.676, CEP. 13.565-905 

February 1, 2008 



Abstract 

We found a new kind of soliton solutions for the 5-parameter family of the potential- free 
Stenflo-Sabatier-Doebner-Goldin nonlinear modifications of the Schrodinger equation. In 
contradistinction to the "usual" solitons like {cosh [P{x — kt)]}~" exp[i{kx — ut)], the new 
Finite-Length Solitons (FLS) are nonanalytical functions with continuous first derivatives, 
which are different from zero only inside some finite regions of space. The simplest one- 
dimensional example is the function which is equal to {cos [j{x — kt)]}^~^^ exp[i(A;x — cot)] 
(with 6 > 0) for \x — kt\ < '/r/27, being identically equal to zero for \x — kt\ > 7r/27. 
The FLS exist even in the case of a weak nonlinearity, whereas the "usual" solitons exist 
provided the nonlinearity parameters surpass some critical values. 
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1 Introduction 



Recently, different authors [||, ||, H, H, ||, ^ ||, ^ |T0|, |TT]] discovered an interesting multi- 
parametric nonlinear homogeneous modification of the Schrodinger equation in the coordinate 
representation. In the most general form this StenfloSabatier-Doebner-Goldin (SSDG) equa- 
tion can be written as (we confine ourselves to the case of a free motion and assume h = m = 1) 

^^ = -1^'^ + (1) 

where the local nonlinear functional is as a linear combination of terms Aip/ip, (Vip/ip)"^, 

jVV'/V'P ^iid their complex conjugated counterparts, so that satisfies the homogeneity 

condition Q{'yilj} = ^liip} for an arbitrary complex constant 7. The specific choices of the 
complex coefficients in the linear combination correspond to the equations describing waves 
in plasmas with sharp boundaries and in nonlinear media [jl], ^ ^, |^. However, trying to 
interpret (|1|) as a quantum mechanical equation one must worry about the conservation of 
probability. For this reason, the functional ^{'ip} was chosen in |^, in an explicit real form: 
fliip} = ^In where T> is the second order differential operator T>f = aA/+b- V/+cV/- V/, 
with real parameters a, b, c. However, it was shown in |1^ that the normalization could be 
saved even in the presence of imaginary (antihermitian) nonlinear corrections of a special kind. 
The most general parametrization was proposed in where ^{ip} was written in terms of 
real and imaginary parts, ^{ip} = R{4'} + ^-^{"^j) ^is follows, 

5 5 

r{^}=dy: a, a, bp] = DY: cjRj m , (3) 

where all the coefficients Xj and cj are real, and the functionals Aj [■?/'] or Rj[ip] are expressed 
in terms of the derivatives of the wave function or in terms of the probability density p = ip^ip 
and the probability current j = (■?/'* V'?/' — ipVip*) /2i: 
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The coefficients Aj and are related as follows, 



Ai = 2c2 
A2 = Ci 
A3 = 2C5 
A4 = C4 

A5 = 2C2 + 2C5 + 



Ci = A2 
C2 = |Ai 

C3 = A5 — Ai — A3 
C4 = A4 

C5 = 4(A5 + A3 — Ai^ 



2^ A ^3 41 

More general homogeneous nonlinear functionals, which include as special cases the nonlocal 



terms proposed by Gisin and by Weinberg [0, were given in |jT^, |T5|, |16|, |T7 



It is not clear, until now, whether nonlinear corrections to the Schrodinger equation of the 
SSDG type have a physical meaning from the point of view of quantum mechanics (possible 
experiments which could verify the existence of such corrections were proposed in |jT8|, |l9l , and 
the relations between the SSDG-equation and the master equation for mixed quantum states 



were studied in [p!5| , [19], ^). Nonetheless, the mathematical structure of the new family of 
nonlinear equations appears rather rich. In particular, studying this family resulted recently in 



discovering the nonlinear gauge transformations |21, 22, 23 



The aim of our article is to show another remarkable property of the SSDG equation, namely, 
the existence of a new type of soliton solutions, which are different from zero in a finite space 
domain even for arbitrarily small nonlinear coefficients. As far as we know, such kind of solitons 
was not discussed earlier. 
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2 Soliton solutions with linear phase 

Looking for a shape invariant solution to the SSDG equation with a linear phase, 

we obtain the following two equations for the real function glx): 



(k-v)^ = D 

9 



l-a)^-ef^V-2;.k.^ = k^(l + r,)-2., 
9 \ 9 J 9 



where the new coefficients are defined as 

a = 2DXi = ADc2, ^ = 2D (A3 + A5) = AD (ca + 2C5) , 

r] = 2D (As - A3 - Ai) = 2Dc3 , /x = 2D (A2 + A4) = 2D (ci + C4) 
A general solution to eq. (13) in the one-dimensional case is 



(4) 



(5) 
(6) 



goix) = {Ci + C2exp[{k -v)x/D]} 



1/2 



with arbitrary constants Ci and C2. However, the function goix) cannot be normalized, thus 
in order to guarantee normalization we impose 

k = v, D = 0, 



i.e. soliton solutions can only exist in the absence of dissipative terms in the Hamiltonian. 
The substitution 

^7(x) = [/(xr, a= 

1 - o- - ^ 

eliminates the nonlinear term (Vg/g^ in eq. (|^), such that 



(7) 



VV-2a^-V/ + 7'/ = 



where 
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2cj-k2(l + r/) 



a 



2/ik 

1-a 



(9) 



Note that 7^ is a free parameter, which may assume both positive and negative values, depend- 
ing on the packet average energy 



/oo 
^*(x,t) 
-00 



gV'(x,t) 
dt 



(ix = uj. 



A general solution to eq. (||) in the one-dimensional case reads 



f{x) = e^^{Cre'^ + C2e-'^), s = ^k^- 7^. (10) 

In particular, 

(a) For 7^ < 0, function ([TOD goes to infinity when x ±00 (if both constants Ci and C2 are 
positive), so a normalizable solution (yf(x) (eq. (|^)) exists only under the condition a < 0, i.e., 
for parameters a and ^ satisfying the inequalities a < 1, ^>1 — aoTa>l, ^<1 — cr(in 
other words, these parameters must be located between the straight lines cr = 1 and a + ,^ = 1 
in the cx^-plane). This means that only strong nonlinearity can give "usual" soliton solutions 
with exponentially decreasing tails, whose simplest representative (/x = 0) reads 

(7,(x) = [cosh(/3x)]-l"'. (11) 

This conclusion agrees with the results of studies [0, |3, , where exponentially confined solitons 
were found for the nonlinear functionals like = aA(ln Similar solutions to the special 

cases of the SSDG equation with complex coefficients were studied in |, A large family of 
exact solutions corresponding to the most general Doebner-Goldin parametrization (@)-(|3) was 
found in 123, E^, p6[. However, that family does not contain the solitons with a linear phase. 



For example, the solution given in has the same amplitude factor as in (|TT]), but its phase 
is proportional to ln[5f(x — kt)], so it does not converge to the plane-wave solution of the linear 
Schrodinger equation when the nonlinear coefficients D and D go to zero. 

(b) For < 7^ < K^, function (^) goes to infinity only for x —>■ +00, while for x —00 it goes 
to zero (or vise versa). In this case, we cannot obtain a normalizable solution in the form (|^) 
for any value of a. 

(c) Quite different situation arises when 7^ > k^, then expression (p!0| ) shall contain trigono- 
metric functions, and (making a shift of the origin, if necessary) we arrive at a solution to eq. 



in the form 



with 7 = y/'-f'^ — n'^ > and 



(12) 



(13) 



1-a-e 

At first glance, we have a problem when / < 0, since function f" is ill-defined in this case 
(unless the exponent a is an integer). But we notice that if a > 1 (i.e. 6 > 0), then the 
function g{x) = [/(x)]° turns into zero together with its derivative g'{x) at 7X = it/2. This 
means that there exists an integrable solution with a continuous first derivative, which is 
localized completely inside a finite domain: 

{cos [^{x-kt)]y^^exp[{l + 5) K{x-kt)+i{kx-ujt)] if \^{x - kt)\ < tt/2 

if \^{x-kt)\>7v/2 

It is remarkable that such a "finite-length soliton" (FLS) exists for an arbitrarily weak nonlin- 
earity, since the requirement 6 > implies the inequalities 



ips{x,t) 



0<^<l-a 



(14) 



which can be satisfied for small values of ^ and a (another possibility is > > 1 — cr, but 
it demands a > 1, meaning a stronger nonlinearity) . In terms of the coefficients Xj and c 
condition (|14[) reads 

D (A3 + As) > 0, 2D (Ai + A3 + A5) < 1 

D (C2 + 2C5) > 0, 8^ (C2 + C5) < 1 . 



It was shown in |jTT[ that the SSDG equation is Galilean invariant provided that (i) ci + C4 = 
and (ii) C3 = 0. In our notation this means /i = 77 = 0. Thus we arrive at the 3-parameter 
family of homogeneous local nonlinear functionals admitting Galilean-invariant and spatially 
confined soliton solutions to eq. ([1|): 
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(15) 
(16) 



Note that parameter v = 2DX2 = 2Dci does not make any influence on the discussed solutions, 
thus only derivatives of the density p and not of the current density j are important for soliton 
solutions. The only crucial parameter is ^, so, the simplest 1-parameter nonlinear functional 
admitting FLS solution reads {u = a = 0), 



Re 



8 V p . 



0<e<l, (17) 



^k7(x,t) 



The explicit form of all FLS-solutions is as follows, 

[/^(x - kt)]^+'' e^C^^-'^"^*) if |x - kt| G 7^(+)(^; 

^^^^ 

if |x-kt| ^7^W(/7) 

where /7(x) is any positive solution to the Helmholtz equation (V^ + 7^) / = with an arbitrary 
real constant 7, and TZ^^\f^) is the internal part of a space region bounded by a closed surface 
(in 3 dimensions) or a closed curve (in 2 dimensions) determined by the equation /-),(x) = (in 
principle, this region may be multi-connected). To avoid any ambiguity, we define the nonlinear 
functional ^{ip} for = assuming VL{'il)}%l) = on such points. Although the solution ([l8| ) 
is non-analytical, it has continuous first derivatives in all points of the space. 

Under the Galilean invariance symmetry (77 = 0), the frequency uj\^^ (eq. (|D) equals 

1, 9 1 9 (1 - ^)^ 

so, the usual dispersion relation of linear Quantum Mechanics (cuk = |k^) is modified by an 
additional constant term (proportional to 7^) that may be interpreted as an "internal energy" 
of the wave packet ([T8|) due to its confinement, whereas k^/2 is the energy of the "center-of- 
mass" motion. For 7 = 0, then V^/o = 0, and considering /o = 1 we obtain a plane wave 
solution to the Schrodinger equation with the wave number k. 

The concrete shapes of the FLS-packets in 2 and 3 dimensions may be quite diverse. The 
most symmetric solutions are given by eq. ([18|) with /^(x) in the form of Jo(7|x|) or sin(7|x|)/|x| 
(in 2 and 3 dimensions, respectively). However, there exist also asymmetric packets with 
functions /7(x) proportional to Jm(7|x|) cosirmp) or j/(7|x|)y/m(''9, v^), where Jm{x) is the Bessel 
function, ji{x) is the spherical Bessel function (proportional to the Bessel function with the semi- 



integral index), and Yim{'&, ^p) is a real- valued analog of the spherical harmonics, ip being the 
usual angular variables. 



3 Discussion 

Although the substitution (|^ was used already in 11 , the existence of the FLS solutions was 
not noticed before, perhaps, because the authors of the cited papers were looking for solutions 
of the stationary Schrodinger equation or for usual exponentially confined solitons. It should 
be noted that the substitution (|^ linearizes only the equation (|^) for the real amplitude of the 
special form of solution (^, but not eq. as a whole. As was shown in [21, 22], the nonlinear 
gauge transformation (NGT) 

ijj' = \ip\exp[i{z*\nijj + zlnijj*)] (20) 

(where z is a complex parameter) transforms any SSDG equation into an equation of the same 
kind, but with another set of coefficients, and all equations can be classified in accordance with 
the possible values of 5 invariants of the NGT family. In the case of eq. (^ with the functional 



Ym the invariants follows (we use the same notation as in Ua 



n = r4 = 0, T2 = -, T3 = -1, 25 = - — . 

8 lo 

If we had = 0, then the whole SSDG equation could be linearized by means of a suitable 
NGT. It is the nonzero value of parameter ^ that prevents the linearization and makes possible 
the existence of the FLS solutions. It is interesting to note in this connection, that various 
special cases of the general SSDG equation were considered before the general structure of the 
equation was found in [§, |ll|, but the coefficients were chosen in such a way that the parameter 



^ was almost always taken equal to zero p7| , |28| , pQ] . The only exception is Ref. ||30|, where the 



only nonzero coefficient in the Aj-parametrization is A5; however, in this case not only ^ 7^ 0, 
but rj ^ 0, too, so the Galilean invariance symmetry is absent. 

Turning to a possible physical meaning of the FLS solutions, we may say that they realize 
the "dream of De Broglie" , in a sense that they permit to identify a quantum particle with a 
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nonspreading wave packet of finite length travelling with a constant velocity in the free space. 
Earlier, the only proposed nonlinear equation that resulted into a nonspreading wave packet 
solution for a free particle was the one proposed by Bialynicki-Birula and Mycielski (BBM), 
with the nonlinear term = —bin lipl"^ (it was shown recently that this term can arise if one 

applies to the SSDG equation a nonlinear gauge transformation with time- dependent coefficients 
p2|). The solitons of the BBM-equation are Gaussian wave packets (gaussons) whose constant 
width is inversly proportional to the nonlinear coefficient b. In contrast to the gaussons, the 
FLS of the SSDG equation have the width 7"^ as a free parameter, independent of the nonlinear 
coefficients. The most attractive feature of FLS solutions is that they exist for an arbitrarily 
weak nonlinearity. Consequently, the superposition principle of quantum mechanics, which is 
verified, from the experimental point of view, with a limited accuracy, does not rule out the 
nonlinear terms like {^n\ ) immediately. On the contrary, new experiments on the verification of 
(non?) linearity of quantum mechanics could be proposed, which would take into account the 
FLS phenomenon. 
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